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Abstract 



Ph . We give a review of the theory of random fields defined on the observable 

r^ I part of the Universe that satisfy the cosmological principle, i.e., invariant 

C^ ' with respect to the 6-dimensional group ^ of the isometrics of the time slice 

of the Friedmann-Lemaitre-Robertson-Walker standard chart. Our new re- 
sults include proof of spectral expansions of scalar and spin weighted ^- 
invariant cosmological fields in open, fiat, and closed cosmological models. 
> 

in 

^ 1 Introduction 

in 

'sf . The space mission Euclid, planned for launch in 2019, is expected, just to mention 

^ i a few, to measure the shapes of 1 .5 billion galaxies. In mathematical terms, Euclid 

will observe a single realisation of a random section of a spin bundle over a ball of 
radius about 10 billion light years. Analysis of Euclid data should base on spectral 
k>i ■ expansions of different cosmological fields. In this paper, we make attempt to 

^ . review the known expansions and to present some new results. 



In Section [21 we review the essential axioms, definitions, and results of the 
Standard Cosmological Model. We consider open, flat, and closed Friedmann- 
Lemaitre-Robertson-Walker cosmological models. The time slice X of each of 
the above models has 6-dimensional group of isometrics $f , where ^ = SOo(l, 3) 
for the open model, ^ = IS0(3) for the flat model, and ^ = S0(4) for the closed 
model. 

Section[3]introduces scalar cosmological fields, both deterministic and stochas- 
tic. In particular, S ub section 13 . 1 [ review s expansions of deterministic cosmological 
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fields in eigenfunction of the ^-invariant Laplace-Beltrami operator. The radial 
coefficients Rkeix) ^^^ RacXx) of the above expansions are used in Subsection l3.2l 
to deduce spectral expansions of ^-invariant stochastic fields. The above expan- 
sions are proved in Theorem [1] similarly for all three models, using tools from 
harmonic analysis. 

In Section |4] we first motivate introduction of spin stochastic cosmological 
fields, using the theory of the effective lensing potential. After short review of 
spin weighted spherical harmonics, we prove Theorem |2] that gives the spectral 
expansion of the convergence, shear, and of the two gravitational flexions cosmo- 
logical fields. 

2 Preliminaries from cosmology 

According to modem cosmology, the Universe is governed by Einstein's equa- 
tions. Solution of the above equations is beyond the capabilities of modem com- 
puters. Cosmologists use another approach. They start from a simple exactly 
solvable model as a background solution and then add perturbations order by or- 
der. 

The background solution is described as follows. The Universe is a 4-dimen- 
sional Lorentzian manifold M satisfying the following two axioms. 

1 . The cosmological principle: at each epoch, the Universe presents the same 
aspect from every point. 



2. The Weyl's postulate after IWevll (Il923h : 



"The particles of the substratum (representing the nebulae) lie, 
in spacetime of the cosmos, on a bundle of geodesies diverging 
from a point in the (finite or infinite) past." 

The "substratum" is the underlying fluid defining the overall kinematics of a 
system of galaxies. The galaxies follow their geodesies. Only one geodesic is 
passing through each point in M, except at the origin. This allows one to define 
the concept of fundamental observer, one for each geodesic. Each of these is 
carrying a standard clock, for which they can synchronise and fix a cosmic time. 

For each fundamental observer, there exists a standard chart {t,x,0,(p) of 
the manifold M centred at the observer in which the metric has the Friedmann- 



Lemaitre-Robertson-Walker fonn: 

diS^ = -c^dt^ + a^{t)[dx^ + fl{x){AQ^ + ^m^eA(p% (1) 

where the function fK depends on the curvature K as 

r(-i^)-i/2sinh(v/=Z;t), K<0, 

fK{x) = lx. K = 0, 

(K-^/^smiVKx), K>0, 

and where t is the cosmic time, x the comoving distance, 6 and cp are spherical 
coordinates on the celestial sphere, c is the speed of light in vacuum, a{t) is the 
scale factor. The comoving distance from the observer to a galaxy is the distance 
the galaxy will have from the observer when it is as old as the observer now. The 
physical distance is the product of the scale factor and the comoving distance. 
We choose the scale factor to be dimensionless and set a(?o) = 1, where to is the 
value of the cosmic time now. Thus, the comoving distance x has the dimension 
of length, while K has the dimension of inverse length squared. The velocity of a 
galaxy with respect to a fundamental observer is called the peculiar velocity. 

The time slice A of the domain D of the above chart is the Riemannian man- 
ifold X with metric a^{t) [dx^ + fxix) (d^^ + sin^ 6 d(p^)]. The geodesic lines of 
fundamental observers are orthogonal to the time slices. The cosmological princi- 
ple implies that each time slice X is the 3-dimensional space of constant curvature. 
When the curvature is negative, the model is called open, the group of isometrics 
of X is ^ = S0o(l,3), and X = lEP, the hyperbolic space. When the curvature is 
equal to 0, the model is called flat, the group of isometrics of X is ^ = IS0(3), 
and X = M^. Finally, when the curvature is positive, the model is called closed, 
the group of isometrics of X is $f = SO (4), and X = S^, the sphere. For all three 
models, the stationary subgroup of the point (?, 0,0,0) is J^ = S0(3). 

The surface S{x) of a comoving sphere of radius x is S{x) = '^7tf^{x)- The 
^-invariant measure on X in the above chart has the form 

dx = f^{x)smOdxd0d(p. 



The comoving distance x runs from to Tt/y/K when ^ > 0, and from to oo 
otherwise. 

The scale factor, a{t), is not observable. However, consider the electromag- 
netic wave of length Ai emitted by a distant galaxy at time ti < to that approaches 



the Earth now. It has another wave length, \q. The ratio Aq/Ai is observable. It is 
denoted by 

Ao 1 , 

T- = 1+z. 
M 

The quantity z is the celebrated redshift. Moreover, 

1 



a{h) 



1+z 



In 1929, Hubble discovered that the redshift is positive for all but few close galax- 
ies. It follows that a{ti) < 1. In other words, the Universe is expanding. The 
galaxies without redshift are those having big enough peculiar velocities to stand 
against the expansion in the standard chart. 
Denote 

Hq = d{to), 

where a dot denote differentiation with respect to t. Hq is called the Hubble con- 
stant. The number 

^' ~ %kG' 
where G is the Newton's gravitational constant, is called the critical density. 

In general, H{z) denote the Hubble constant as measured by an imaginary 
astronomer working at redshift z- The expansion of the Universe is governed by 
the Friedmann equation: 

where the parameters are as follows. 

SIm is the density parameter for non-relativistic matter (which moves with 
speed much less than c with respect to the fundamental observer). We have 

i^M = , 

Pc 

where p^jo is the density of the non-relativistic matter now, att — to. 
Q./^ is the vacuum energy parameter. 



Ac 



j^ 2 



'0 

where A is the cosmological constant. 



Parameter 


Value 


Ho 


67.80±0.77km/s/Mpc 
Q 315+0.016 

W.-3iJ_QQJg 

-0.0010+Sg 



Table 1: Cosmological parameters 



Qic is the curvature parameter. 



Qk^ 






Finally, Q.r is the radiation density parameter. We have 



where pRtQ is the density of the radiation now. When z = 0, we obtain 



The values for the above parameters, adapted from Efstathioul (120131) . are 
shown in Table [B where 1 Mpc = 3.09 x 10^^ m, and the quoted errors show 
the 68% confidence level. The contribution from radiation is negligibly small. 

The comoving distance % from the fundamental observer at z = to an object 
with redshift z is 



Z(z) 



du 



Ho Jo y/aR{i+u)^+aM{i+u)^+aK{i + u)^+aA 



The look-back time ti to the above object is the difference between to and the cos- 
mic time at the moment when photons were emitted by the object. It is calculated 
as 



tL 



1 



dw 



Ho Jo [i + u)^Q.R{l + uf + Q.M{l+uy + Q.K{\+uY + Q.K 



In Friedmann equation, a = is reached after a finite time. Near a = 0, the 
Hubble parameter //(z) explodes. This point is therefore called Big Bang. 



3 Scalar fields 

3.1 Deterministic fields 

Let / : M — 7- M be a deterministic scalar cosmological field. On each time slice X, 
it is natural to decompose the restriction of / to X in eigenfunctions of ^-invariant 
differential operators. 

Let d{x,y) be a distance between the points x and y EX. In all three models, 
X is a two-point homogeneous space. That is: for any two-point pairs x,y EX, 
x' ,y' eX satisfying <j(x,v) —d{x',y '), there exists an isometry g G ^ such that 



gx = x', gy — y'. By (lHelgasonl . 1 19841 . Chapter 2, Proposition 4.1 1), all ^-invariant 
differential operators are polynomials in the Laplace-Beltrami operator A^: of the 
Riemannian manifold X. In what follows we call A^: just the Laplacian. We have 

1 d / ., , d \ 1 



fl{x)dxy''''''dz) fl{x) 

where 

1 d f 2 d 



A=^^ 0^^ + 



is the Laplacian on the sphere §^. 

The eigenvalue problem for the Laplacian A^: is traditionally written as the 
Helmholtz equation 

AKQi,{x) = -{k^~K)Qi,{x), 

where k is called the wave vector and labels the modes, Q]i{x), and where k is 
called the wave number. The wave numbers take the values in [0, °°) ioxK<0 and 
in the set { ( ft) + 1 ) \fK : ft) G Z+ } otherwise, where Z+ is the set of all nonnegative 
integers. Different authors use different normalisation of modes. 

The Helmholtz equation may b e solved by the metho d of separation of vari- 
ables. We use the normalisation by IPeter and UzanI (|2009|) : 



QkUX. e, (?) = {2nf/^RM{x)Yim{e, (p), 
where Y£fn{d, (p) are spherical harmonics, the eigenfunctions of A: 

AYe^ = -£{£+\)Ye^, £>0, -£<m<i. 
Spherical harmonics are related to associated Legendre polynomials Pf through 

Ym{9,cp) = J 4^(^ + ^), PAcose)e r 



The radial part of the Helmholtz equation becomes 



__l d_ 



flil) 






+ 



k^-K 



+ i: 



m) 



Rm{x)=0. (2) 



When i^ < 0, put ft) = kj^f^-K and 






The solution to Q has the form 



Ru{x) 






— 1/2— £ 

where ^_i/2j-ico is the associated Legendre function of the first kind. 
When ^ = 0, the solution to ^ has the form 



(3) 



where ji is the spherical Bessel function. The dimensionless distance in this case 
is r = kx. 

Finally, when if > 0, put (0 = k/\fK- 1 and 



n=0 

Note that M^^i = for £ > CO. The solution to Q has the form 



Rcoeix) 



nMai 



2(£0 + l)2sin(v^Z) ^/^ 



pJJ^ZUcos{VKx)), 



CO 



(4) 



see (lAbbott and Schaefeilll986[ Equation (A21)). 

For K <0, the expansion of a function / G L^(X, dx) into modes has the form 

/(Z,e,(p) = (2;r)-3/2£ £ / k^ff,„{k)Qke,n{X.e,cp)dk. (5) 



For ^ > 0, the above expansion reads 

(T^ \ 3/2 oa (O l 

^ EEL {o)+\ffumkUx.Q.<(>). 

where the Fourier coefficients, fimik), are given by 

fUk) = {2n)-^/^ [ fiX.e,cp)QMmiX.O,(p)fiix)smedxded(p. 



IX 

Assume that ^ < and / depends only on X- Then we have 



1/2+ico^ 

sin(V^Z^) 



p-;e,. (coshr)--/ ' ^^"^^^^ 



Using the formula 



-l/2+ico^ Y ;rsinhr G) 

we obtain 



and dS]) becomes 



fcsmh(V=rz) 



/(Z)^y^r/oo(^j^^^^^d*, (6a) 

/v^TJo ^smh(v-^Z) 



where 



/»(*)= 'r/(z)-^4^5te)dz. (6b) 

2y/7t Jo fcsmh(v-^Z) 



Similar calculations show that in the case of ^ = we have 

1 Tr n.^MkX),2 



where 



nV2 Jo kx 



Mk) = -^ r f{xf^^S{x)dX, (7b) 

ny/2 Jo kx 



while in the case of A" > we have 



/U) = ^w^ E /oo(cQ) 7 ; a to+1)^ 

2v^c^o w + l sinV^;t 



where 



/oo(w) 



2^/% Jq 






(8a) 



(8b) 



The above expansions, with various normalisations, were cons i dered bvlFock 



(ll935h.lSchrodinged(ll939L 1957h . Lifshitz and Khalatnikovl(ll963h . lHarrisonl(ll967h . 
Abbott and Schaefen (|l986h . among others. 



3.2 Stochastic fields 

Let /: Mxfl— j-Cbea stochastic scalar cosmological field defined on a prob- 
ability space (Q,5^, P). We assume that / has a finite variance, E[|/(f,jc)p] < oo 
and that / is mean-square continuous, i.e., the map M — )■ L?-{Q.), {t,x) i— )■ f{t,x) is 
continuous, where L^{Q) is the Hilbert space of all random variables with finite 
variance. 

Consider the restriction f{x) of the field f{t,x) to the time slice X. By the 
cosmological principle, the field f{x) must be ^-invariant, i.e., for any positive 
integer n, for any distinct points x\, . . .,XnEX, and for any gE'i^, random vectors 
(/(xi), . . .,f{xn)Y and {f{gx\), . . .,f{gXn)Y niust have the same distribution. 

In particular, the expected value E[/(a:)] is a constant, say a G C, and the 
autocorrelation function 



R{x,y) = £[{f{x)-a){f{y)-a)] 

is ^-invariant, i.e., R(gx,gy) = R{x,y). Such fields are called wide-sense homo- 
geneous. Note ihdXR{x,y) is a positive-definite function. 

The spectral theory of wide-sense homogeneous randor n fields on topo logi- 
cal groups and homogeneous spaces has been developed by lYaglomI (1196 lb . He 
proved spectral expansions of wide-sense homogeneous random fields on sepa- 
rable topological groups of type I and their homogeneous spaces. A topological 
group of type I is such a group ^ that every un itary representat ion of ^ generates 
an operator algebra of type I (see, for example, Naimarkl (| 19721) ). All three groups 
S0o(l,3), IS0(3), and S0(4), are of type I. 

Moreover, the above groups share the following property. Let U be an irre- 
ducible unitary representation of the group Sf . Consider the restriction of U to 



the compact subgroup J^ = S0(3) . The above restriction is equivalent to a direct 
sum of irreducible unitary representations of J(f^. The multiplicity of the trivial 
representation of the group ^ in the above sum is equal either to or to 1 . 

Let ^^ be the set of all equivalence classes of irreducible unitary representa- 
tions of the group ^ for which the above multiplicity is equal to 1. Let U^ G ^^, 
and let Hqq be the one-dimensional complex Hilbert space in which the trivial rep- 
resentation of the group J^ is realised. Let e be a vector of unit length in Hqq. It 
is easy to see that the function 

does not depend on the choice of the vector e. It is called the zonal spherical 
function of the group ^. 

Let n: ^ -^ X be the natur al projection, let x, y e X, and let ^i G n^^{x). 



gi&Tt (v). By (lYaglomlll96lL Theorem 6'). formula 



R{x,y) 



^cofeVOdvlw) 



determines a one-to-one correspondence between the autocorrelation functions 
R{x,y) of wide-sense homogeneous random fields on X and finite measures v on 
^jf . When X is a two-point homogeneous space, the zonal spherical function 
^m{S2^Si) depends only on r, the distance between x and y. 

The set ^^^ and the zonal spherical fu nction for the case of Z = H were de- 
termined by iGel'fand and NaimarkI (|l946|) . There ex ist different par ametrisations 
of the above set. We will use the parametrisation by iKostantI (Il975|) . The set ^j^ 
includes the principal series parameterised by ft) G [0, -foo), the supplementary 
series parameterised by OJ G (0, i), and trivial representation with (o = '\. We have 



^coW 



sin(ft)r) 
ft)sinhr 



and 



R{x,y) 



sin(ojr) 



dv(ft)) 



(9) 



ft)sinhr 

where the in tegral is taken over the set (0, i] U [0, °°) . Equation Q was obtained by 
KreinI (119491) for the case of a finite-dimensional hyperbolic space H" . To compare 
^ with (|6al) . we choose r in such a way that or = kx, or r = \^—Kx- We obtain 



Rix) 



1 



y/^Ksm{kx) 



dv{k). 



10 



Compare the last display with (|6al ). Apart from normalisation, we note the 
following difference. The expansion of a square-integrable function includes only 
the representations of the principal series of the group Sf (in physical terms, sub- 
curvature modes, when the wavelength is less than the radius of curvature). On the 
other hand, the expansion of a positive-definite function includes the representa- 
tions of both the principal and the supplementary series (sub- and super-curvature 
modes, when the wav elength is greater than the ra di us of curvat ure). This dif- 



Lvthl (119951) describes this 



ference was noted by iLyth and Woszczynal (|1995h . 
situation as follows. 



Mathematicians have known for almost half a century that all modes 
must be included to generate the most general homogeneous Gaus- 
sian random field, despite the fact that any square integrable /Mncrion 
can be generated using only the sub-curvature modes. The former 
mathematical object, not the latter, is the relevant one for physical 
applications. 

The set '^j^ and the zo nal spherical functions for the case of X = 
determined by .Krein (|l949|) . We have '^jf = [0, o°) , 



were 



4) 



sin( cor) 



a\ 



cor 



To compare this with (TTal ). choose CO = k and r = X- We obtain 



Rix) 



kX 



■dv{k). 



CartanI (|1929|) determined the set ^^ and the zonal spherical functions for the 
case of X = S-^. In this case ^j^ = Z_|_, 



To compare this with (l8al ). put r = 

oo 

R{X) = E v(co 



sin(((0 + l)r) 



C0 + l)sinr 
Kx and recall that 0) = k/y/K- 1 . We obtain 
sm{{co+l)VKx) 



m^o (w + l)sin(v^^) 

where v(ft)) > and Em v(ft)) < 0°. 
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Equations Q, ©, and ([8]) have the following group-theoretical interpretation. 
Consider the set of double cosets J(f\^ /J^f. When K <0, this set is the inter- 
val [0,00), otherwise this set is \(),K\/K]. The measure S{x)dx is a unique up 
to a constant multiplier ^-bi-invariant measure on J^\S j J^ . For any choice 
of the above multiplier, say C, there exists a unique measure on Sfj^ called the 
Plancherel measure. The Plancherel measure is characterised by the following 
property: equations ©-([S]) determine an isometric isomorphism between the 
space of J^-bi-invariant square-integrable functions on 5^ with measure CS{x) dx 
and the space of square-integrable functions on 5^^ with Plancherel measure. 
When ^ < 0, the measure k dk is proportional to the Plancherel measure on ^^. 
Note that the Plancherel measure is equal to zero for the supplementary series 
0) G (0,i) and the trivial representation CO = i. When ^ > 0, the Plancherel mea- 
sure of a point CO E '^x is proportional to (ft) + 1)^. Equations (|6bl) . (iTbl) . and (l8bl) 
define the spherical Fourier transform, while equations (l6al) . (TTal) . and (|8al) define 
the inverse transform. 

To find the spectral expansion of the random field /, we have to find a mea- 
surable space (y, QJ, \x) and a function /? : X x V — )• C such that 

R{x,y) = [ h{x,v)hMd^{v). (10) 

Jv 

In this case there exists a complex- valued scattered random measure Z on V with 
/i as its control measure, i.e., for any A, 5 G QJ we have 



E[Z{A)Z{B)]=^{AnB) 
such that 



f{x)= [h{x,v)dZ{v) 
Jv 



Note that i{h{x, v) satisfies (fTOl) . then h{x, v)e^'f^^'' also satisfies the above equation. 
For any CO G ^jf^, consider the restriction of the representation [7® to ^. In 
the cases of the open model with co ^ i and the flat model with ft) 7^ 0, the above 
restriction is equivalent to the direct sum of the irreducible unitary representations 
r^ of the group J^ for all t G Z-|- acting in Hilbert spaces H^. In the case of the 
closed model, the above restriction is equivalent to the direct sum of the represen- 
tations D^ , < £ < ft). The restriction of any representation D^ to the subgroup 
S0(2) C J^ is equivalent to the direct sum of the representations (p 1— )■ e™*'' of the 
group SO (2), ~£ <m< £. We represented the space Hco of the representation [7® 



12 



as the direct sum of one-dimensional subspaces //^„. Choose a vector e^^ G i/^^ 
of unit length in such a way, that the entries of matrices 

are Wigner D- functions: 



Dl((p,0,y)=.-H-^+"V^)(-ir /^^^+^|;^ 



+ n\ f in \^ .^£-s+n^^.2s-m-n, 



mm{£+m,£+n} 

X y |-'"j| - " )(-l)^--'+«cot^-'-"'-''(0/2). 

m , 1 V ■^ / \s — m — nl 

(11) 

Let g = kirk2 be the standard polar decomposition of the elements of the 
group ^. Choose gi = kiVi G K^^{xi), i= 1,2. Then we have 






By matrix multiplication 

E 






The matrix entry tZ/^^-ool'"!) i^ called the associated spherical function. For all 
three models, they are different from only when m = 0. Therefore, 



Rix,y) = '£DUk2'ki) t/,Vo(n)t/5,;o,o('-2)dv(«). 



Again by matrix multiplication 



e m=-£ -^^JT 

Wigner D-functions are related to spherical harmonics by 



2£+l 
4;r 



Y£mid,(p) = \l-^D'y^Q{(p,e,0). (12) 
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Therefore we have 



1 ^ 



X I, f^w,o(n)C/,Vo('-2)dv(co), 
where 

}^lu K = 0, 

and (;fi, 0i, 91) (resp. (xi^ 02, 92)) are coordinates of the point x (resp. y) in the 
standard chart. 

In the case of the open model we have 

V = {{{0,i]U[0, 00) J, m):ieZ+,-e<m<i}, 

i.e., the union of countably many copies of the set (0, i] U [0, 0°) indexed by pairs 
(£, m), 5J is the a-field of Borel sets in V, and the restriction of the measure jU to 
any copy is equal to v. 



Th e associated spherical functions for this case were calculated by IVilenkin 
(Il958h : 






Substitute this formula to (fT3] ). We obtain 

00 £ 

R{x,y) = 2n^Y. E >£m(0i,<Pi)Jfm(02,<P2) 

f=Om=-£ 



ir(i«)|2p_;;2VJ'^o^Mv/=^zi))^_i;2Vco('^°^h(v/=Zzi)) 

X / , dv(£o) 

\r{ico-i)\\ smh{V^Kxi)smh{y/^Kx2) 



Using the formula r(z) = (z — l)r(z — 1) and ^, we have 

E 



•x> i poo 

fiX,d,(p) = 2y^Y. E YUe,(p) RMix)dZMik). 

ij—r\ , — ij Ji\/—K 
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The case of the flat model is similar, but this time 

V = {{[0,oo)J,m):£eZ+,-i<m<£}. 
The associated sphe rical functions for this case were known in implicit form since 



Erdelvi et al. 



9531) . The group-theoretical interpretation has been proposed by 



IVilenkini (Il957h . We have 

f^W,oW=i'v/2ZTT7>(a)r), 
and the spectral expansion has the form 

Finally, in the case of the closed model 

y = {(£t),£,m): ft)GZ+,0<£<a),-£<m<£}, 

2J is the a-field of all subsets of V, and iJ.((0 , i,m) = v{(o) . The associated spher- 
ical functions are known since Erdelvi et al. (1953) and have the form 



{(0~i)\{2i+l] 



where C^^^£(cosr) are Gegenbauer polynomials. Using the formula 

C^(cosr)- ^^^ + ^^~^^'^^"'^'"'V /^-^ 

the value ft) = k/y/K— 1, and dU, we obtain 

U^O;0,oir) = VliTlRkiiX) 
and 

CO 0) £ 

(0=0 £=0 m=-e 
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In what follows we suppose that the measure V is absolutely continuous with 
respect to the Plancherel measure. In other words, 



Rix) 






where P{k) is the density of the measure v with respect to the Plancherel measure. 
P{k) is called the spectral density by mathematicians and the power spectrum by 
cosmologists. The power spectrum P{k) is the spherical Fourier transforms of 
the autocorrelation function R{x), while the latter is the inverse spherical Fourier 
transform of the power spectrum. We have the following theorem. 

Theorem 1. The spectral expansion of the '^ -invariant random field f takes the 
following form. In the case of the open model 

/U,0,(p) = 2v^£ £ i£m(0,9)/ RkMkVmdWUk). (14) 
In the case of the flat model 

£=Om=-e -^0 

Finally, in the case of the closed model 



co=0e=0 m=-£ 

where Wi^ is a sequence of uncorrelated identically distributed complex-valued 
scattered random measures on [0,°°) with Lebesgue measure as their common 
control measure, and where W^Qfj^ is a sequence of uncorrelated identically dis- 
tributed random variables with zero mean and unit variance. 



4 Spin fields 

In order to motivate the introducing of spin fields, consider the following example. 
Let p{t,x), X = (;^, 0, (p) G X be the density of the non-relativistic matter. Cur- 
rently, only the baryonic part of p is observable. Define the fractional overdensity 
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as 

pM,t 

The Newton's gravitational potential, ^, is the solution to Poisson's equation 

AK^{t,x) = — ^p(?,x). 

Consider a source of electromagnetic waves at the comoving distance X with an- 
gular coordinates 6 and (p. The ejfective lensing potential of the source is defined 
as 

yriX.e.cp) = -^ r ■^''[^~'U {to-tL.r,0,cp)dr. 

jK[X) Jo fKin 

This quantity is not observable. However, the convergence K = (l/2)9*9i//' and 
the shear 7= (1/2)5^1// are observable now, while the first gravitational flexion 
F = (1/2) 59* 5i// and the second gravitational flexion G = ( 1 /2) 5^ i/A may become 
observable in the near future. 

The symbol 5 denote the spin raising operator. There exist several equivalent 
definitions of this object. Consider the one connected to representation theory. 

Let 5 G Z be an integer number. The map 5 1— )> ((jO 1— )> e"'*'), (p E S0(2) de- 
termines a one-to-one correspondence between Z and the set of the irreducible 
unitary representations of the group Jif = S0(2). Consider the Cartesian product 
S0(3) X C and introduce the following equivalence relation: the points (^i,zi) 
and (^2722) are equivalent if there exists an element h = hq, E Jif such that g2 — 
g\h and zi = e^"*''zi. Call the quotient space E^. The projection map 7r{g,z) ~ 
g,^ maps Es onto the homogeneous space SO(3)/SO(2), which is the two-di- 
mensional sphere, S^. The inverse image of each point in S^ is a copy of C. The 
triple {Es, %, §^) is a homogeneous line bundle over §^. 

Let / : §^ — 7> Es be a section of the above bundle, i.e., ;r o / is the identity map 
on S^. The group SO (3) acts on Es by g{gQ,z) = (ggo.z). This action identifies 
the fibers over any two points of the base space. Therefore, we can define the 
square-integrable sections by the condition 

i'2k pK 

/ / |/(0,(p)|2sineded(p<oo. 
JO Jo 

The elements of this space are called the square-integrable functions of spin s. 

The induced representation U of the group S0(3) acts in the space L2{Es) of 
the square-integrable sections by 

U{g)nn) = f{g-'n), gE S0{3), nES\ 
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The Frobenius reciprocity states that the multiplicity of the irreducible representa- 
tion U(_ in U is equal to the multiplicity of the representation sinU^. The latter is 
equal to 1 if ji'l < £ and to otherwise. It follows that the induced representation 

U is equivalent to the direct sum of irreducible components t/i^, £ = |5|, 1^1 + 1, 

The basis in the space L2{Eq) consists of spherical harmonics (fT2] ) (of spin 0). 
The basis in Li^Es) consists of spherical harmonics of spin s: 



The last display and (fTTI) imply the following formula: 



JUe,cp)-e (-1) W 4n{i + s)Ki-s)l ''" ^^^^^ 

voin{£+m,£-s} /« _ \ / « , \ 

X E 1 )( )(-l)^-"-^cot2"-™+-'(0/2). 

■7^ .\uj \u — m + sj 

M=max|0,m— ij v / v / 

Consider the operator 5 (in fact, a family of operators) acting onto a spherical 
harmonic of spin s by 



This action may be extended by linearity to the linear subset ^g of the space 
L2{Es) that consists of the functions 5/, for which their Fourier coefficients 

sfim= I sfiO, (p)sYe^{e,(p) sin Ode d(p 
satisfy 

00 

£ {£-s){i + s+l)l.fM\^<oo. 

i=\s\ + l 

It is easy to check, that 

r. d i d 

a = scotO 



do smO d(p' 
The conjugate operator, 5*, acts onto a spherical harmonic of spin s by 
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By obvious reasons, 9 is called the spin raising operator, while 5* is called the 
spin lowering operator. Moreover, the restrictions of these operators to the space 
H(, where the irreducible representation U^ acts, are intertwining operators be- 
tween equivalent representations. 

Assume that the effective lensing potential of the sources of radiation is a 
single realisation of a random field. By the cosmological principle, the above 
field must be ^-invariant. In particular, its spectral expansion is determined by 
Theorem[T] Consider the case of the convergence field in an open model. Assume 
that the power spectrum P{k) satisfies the following condition: 

££(£+l)(2£+l) / \RM{x)\'k^P{k)Ak<^. 
f=i Jo 

Apply operator (1/2)5*3 to the expansion (fl?)) term by term: 

I 

E 



D—\ ^ — _ll Jo 



The series in the right hand side converges in mean square and determines the 
spectral expansion of a random field. Note that the above expansion cannot 
be written in the form of ([14]) . Therefore, the random field K{x,6,(p) is not 
S0o( 1 , 3) -invariant. However, the restriction of the field K to the sphere of a fixed 
comoving radius Xo centred at the fundamental observer is an isotropic random 
field on the sphere. In other words, the above restriction is an isotropic random 
section of the line bundle (Eq, ;r, S^). 

The heuristic arguments of the previous paragraph suggest the following. The 
restriction of the convergence field to the sphere of a fixed comoving radius Xo 
centred at the fundamental observer is an isotropic random section of the line 
bundle {Eo,n,S^), that of the first gravitational flexion field is an isotropic ran- 
dom section of the line bundle {Ei,n,S^), that of the shear field is an isotropic 
random section of the line bundle (£'2, 7t,S^). Finally, the restriction of the second 
gravitational flexion field to the sphere of a fixed comoving radius Xo centred at the 
fundamental observer is an isotropic random section of the line bundle (£3, ;r, §^). 
This means that the expected value E[X(;fo5n)], X E {fc, 7,F, G}, n G S^, is con- 
stant, while the covariance function E[X {Xo,ni)X (xo^'^i)] satisfies the following 
condition: 

E[X{xo,gni)X{xo,gn2)]=E[X{xo,ni)X{xo,n2)], geS0{3). 
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Spectral expansion of an i sotropic random field in the trivial bundle (Eq, ;r, S^) 
goes back to Obukh ovl(ll947b. that in the bundle (£'2, ;r, S^) is known to cosmol- 
ogists since Zaldarriaga and Seljakl(|l997|) . The rigourous r nathematical theory of 



isotrop ic r andom fields in vec tor bundles was proposed by iGeller and Marinucci 



op 

(|2010|) and lMalyarenkol (|2011t) . In particular, the spectral expansion of an isotropic 
random section Z(0, (p) of the line bundle (Es, TT, §^), s eZ, has the form 



where E[aj^£OT] = unless s = £ = and 

00 

£(2£+l)Crf<oo. 

£=s 

It follows that there exists a sequence of uncorrected mean square continuous 
stochastic processes as£mix) such that 



£=s m=—£ 

Moreover, E[as/,n{x)] = unless s = i = m = and 



with 



^[as,£miXl)asI'm'{X2)] = %5„m'QKZl, Z2) 



£(2£+l)Q(X,Z)<-. 



(15) 



(16) 



(17) 



As usual, the comoving distance x runs from to n/y/K when ^ > 0, and from 
to oo otherwise. 

The correlation function of the random field (flST) is 



R{Xh0l,(PhX2,d2,(p2) = ^[X{Xh0l,(Pl)X{X2,d2,(p2)] 



£Crf(Zl,Z2) £ sY£„r{dh(Pl)sY£Ud2,(p2)- 

£=s m=—£ 
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To simplify this formula, use the addition theorem for spin weighted spherical 
harmonics from/ iDurrerl (|2008[ equation (A4.57)): 



An 



2^q7Y ^ sYM'{d2,(p2)-mYM{0i,(Pi)=sYe,m{l3,a)e ''^, 



where (a,/3, 7) are the Euler angles of the rotation which first rotates (62^ (pi) into 
the north pole and then north pole into (0i, (pi). In fact, /3 is the angle between 
(01, (pi) and (027 92)5 Of is the angle between the great arc connecting the above 
two points and the meridian passing through (Oi.cpi), and 7 is the angle between 
the great arc connec ting the above two points and the meridian passing through 
(02, 92)7 see Fig. 1 in iNg andLiul(|1999() . Substitute the value m = —s, use formula 



,F,„(i8,a) = e-""V,„(i8,0) 



(18) 



and interchange indices 1 and 2. We obtain 



2^ sYlm ( 01 , 9l ) sYhn ( 02 , 92 ) 



2i+l 
Atz 



sYe-s{^.0)e 



-is(a+y) 



and 



^(zi,0i,<Pi,;t:2,02,<?>2) = :r^£Q(zi,;k:2)\/2?TT.y£-.(j8,o)e 



-is{a+y) 



(19) 

To calculate Csi{X\,X2), put 0i = (P2 = and 02 = j8 in the last display. In 
other words, let (0i, ^i) be the north pole. Then a = 7 = and 

-i 00 

^(Zi,0,0,Z2,i8,0) = -^£c,,(;t:i,Z2)v/2£TT..y/-.,(/3,0). 

The spin weighted spherical harmonics sYimiQ^^, ^> 1*1, —i<m<£ form an 
orthonormal basis in the space of square integrable functions on the sphere §^. It 
follows from (fTSi) that the functions y/lTrsYi-si^i^) form an orthonormal basis 
in the space L^([O,;r],sin0(i0). Multiply both hand sides of the last display by 
sY(^sil^:0) and integrate over [0, ;r] with respect to the measure sin/3dj8. We 
obtain 



CsiiXhXi) 



A%^l^ 

VilTTJo 



^(Zi,0,0,Z2,/3,0),y,_,(/3,0)sin/3d/3. (20) 
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Multiply both hand sides of equation (fT5] ) by sYf,n{6,(p) and integrate over §^. 
We obtain 

as/mU)= I X{x,e,(p)sYem{0, (p) sin0 d0 d(p. 

It follows from mean square continuity that 

limE[|a,,^^U)P]=0, (^,£,m)^ (0,0,0), 

therefore 

C,,(0,Z)=0, i^O. (21) 

We arrive at the following theorem. 

Theorem 2. The spin cosmological fields have the form (fT5l) . where 5 = /or 
?/?e convergence, s = I for the first gravitational fiexion, s = 2 for the shear, and 
s = 3 for the second gravitational fiexion. The sequence of stochastic processes 
as/mix) satisfies (fT6l) and (fTT]) . The correlation function of a spin cosmological 
field has the form (fT9l ). where CsfiXi^Xi) is a sequence of nonnegative-definite 
kernels satisfying (flTI) and (|2T]). T/ze kernels CseiXi^Xi) ^<^y be calculated by 
(EQl). 



Theorem|2]in the case of i' = is due to I Yadrenkd (| 1 9 8 3|) . 



References 

L. Abbott and R. Schaefer. A general, gauge-invariant analysis of the cosmic 
microwave anisotropy. Astroph. J., 308:546-562, 1986. 

E. Cartan. Sur la determination d'un systeme orthogonal complet dans un espace 
de Riemann symmetriqueclos. Rend. Circ. Mat. Palermo, 53:217-252, 1929. 

R. Durrer. The cosmic microwave background. Cambridge University Press, 2008. 

G. Efstathiou. Planck 2013 results. XVI. Cosmological parameters. Astronomy 
and Astrophysics, 2013. arXivl303.5076vl [astro-ph.CO]. 

A. Erdelyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi. Higher transcenden- 
tal functions, volume II. McGraw-Hill, New York, 1953. 

V. Fock. Zur Theorie des Wasserstoffatoms. Zeitschrift fUr Physik, 98:145-154, 
March 1935. doi: 10. 1007/BFO 1336904. 

22 



I. M. Gel'fand and M. A. Naimark. Unitary representations of the Lorentz group. 
Acad. Sci. USSR. J. Phys., 10:93-94, 1946. 

D. Geller and D. Marinucci. Spin wavelets on the sphere. J. Fourier Anal. AppL, 
16:840-884, 2010. 

E. Harrison. Normal modes of vibrations of the Universe. Rev. Mod. Phys., 39(4): 
862-882, 1967. 

S. Helgason. Groups and geometric analysis. Integral geometry, invariant dif- 
ferential operators, and spherical functions. Pure and applied mathematics. 
Academic Press, Orlando, 1984. 

B. Kostant. On the existence and irreducibility of certain series of representations. 
In I. Gel'fand, editor. Lie groups and their representations, pages 231-329. 
Halsted Press, New York, 1975. 

M. Krein. Hermitian positive kernels on homogeneous spaces. Ukrain.Mat. Zur- 
nal, l(4):64-98, 1949. 

E. Lifshitz and I. Khalatnikov. Investigations in relativistic cosmology. Advances 
in Phys., \2:\S5-249, 1963. 

D. Lyth. The Grishchuk-Zeldovich effect in the open Universe. Annals N. Y. Acad. 
Sci., 7 59:101,1995. 

D. Lyth and A. Woszczyna. Large scale perturbations in the open Universe. Phys. 
Rev D, 52:3338-3357, 1995. 

A. Malyarenko. Invariant random fields in vector bundles and application to cos- 
mology. Ann. Inst. H. Poincare Probab. Stat., 47(4): 1068-1095, 2011. doi: 
10.1214/10- AIHP409. 

M. A. Naimark. Normed algebras. Wolters-Noordhoff Ser. Monogr. Text. Pure 
Appl. Math. Wolters-Noordhoff Publishing, Groningen, third edition, 1972. 

K.-W. Ng and G.-C. Liu. Correlation functions of CMB anisotropy and polariza- 
tion. Int. J. Mod. Phys. D, 8(l):61-83, 1999. 

A. M. Obukhov. Statistically homogeneous random fields on a sphere. Uspehi 
Mat. Nauk, 2(2): 196-198, 1947. Russian. 



23 



p. Peter and J.-P. Uzan. Primordial cosmology. Oxford University Press, Oxford, 
2009. 

E. Schrodinger. The proper vibrations of the expanding Universe. Physica, 6: 
899-912, 1939. 

E. Schrodinger. Expanding Universes. Cambridge University Press, Cambridge, 
1957. 

N. Ya. Vilenkin. The matrix elements of irreducible unitary representations of the 
group of real orthogonal matrices and group of Euclidean {n — 1) -dimensional 
space motions. Dokl. Akad. Nauk SSSR (N.S.), 113:16-19, 1957. Russian. 

N. Ya. Vilenkin. The matrix elements of irreducible unitary representations of a 
group of Lobachevsky space motions and the generalized Fock-Mehler trans- 
formations. Dokl. Akad. Nauk SSSR (N.S.), 118:219-222, 1958. Russian. 

H. Weyl. Zur allgemeinen Relativitatstheorie. Physikalische Zeitschrift, 29:230- 
232, 1923. 

M. I. Yadrenko. Spectral theory of random fields. Translat. Ser. Math. Eng. Opti- 
mization Software, Publications Division, New York, 1983. 

A. M. Yaglom. Second-order homogeneous random fields. In Proc. 4th Berkeley 
Sympos. Math. Statist, and Probab., volume II, pages 593-622, Berkeley, CA, 
1961. Univ. California Press. 

M. Zaldarriaga and U. Seljak. An all-sky analysis of polarisation in the microwave 
background. Phys. Rev. D, 55(4): 1830-1840, 1997. 



24 



